Dynamics of coherent polaritons in double-well systems 
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We investigate the physics of coherent polaritons in a double-well configuration under a resonant 
pumping. For a continuous wave pump, bistability and self-pulsing regimes are identified as a 
function of the pump energy and intensity. The response to an additional probe pulse is characterized 
in the different cases and related to the Bogoliubov modes around the stationary state. Under a 
pulsed pump, a crossover from Josephson-like oscillations to self-trapping is predicted for increasing 
pump intensity. The accurateness of the effective two-mode model is assessed by comparing its 
predictions to a full solution of the non-equilibrium Gross-Pitaevskii equation. 

PACS numbers: 71.36.+C, 71.35.Lk, 42.65.-k, 03.75.Lm 



I. INTRODUCTION 

Semiconductor microcavitics in the strong coupling 
regime are particularly well suited to study the physics 
of dilute Bose gases in a solid state context^* 2 -! 3 -! 4 -. The 
elementary excitations of the system consist of polari- 
tons, i.e. a superposition of a cavity photon and an exci- 
ton which at low excitation levels satisfy Bose statistics. 
Their photonic component guarantees that a large degree 
of spatial coherence is maintained in spite of disorder ef- 
fects, while the excitonic one provides strong mutual in- 
teractions. Differently from most other examples of Bose 
gases such as liquid 4 He and ultracold atoms, a polari- 
ton gas is an intrinsecally non-equilibrium system, whose 
properties can dramatically differ from the corresponding 
ones of systems at thermodynamical equilibrium^! 7 ^. 

Recent advances in the semiconductor fabrication tech- 
nology, have made it now possible to design polariton 
traps with a high flexibility in the shape and the depth 
of the trapping potential ^ 10 ' 11 : 12 ' 13 : 14 ' 15 . From this per- 
spective, double-well potentials show a particular in- 
terest as they provide a way of investigating the well- 
known Joscphson effec t 16 : 17 ' 18 in completely new non- 
equilibrium regimes. Some preliminary work for the case 
of non-resonantly pumped polariton condensates has re- 
cently appeared in 8 -, while many authors have consid- 
ered similar effects in a variety of different optical sys- 

tem9 19 1 20 i 21^ 

In the present paper, we will concentrate on the case of 
resonantly and coherently pumped double-well polariton 
traps obtained by lateral patterning of a planar microcav- 
ity as experimentally done i n 12 ' 13 : such a configuration 
allows not only for selective addressing and diagnostics 
of the two spatial modes, but also for a relatively easy 
time-resolution of the Joscphson dynamics on a picosec- 
ond scale. The mean-field calculations of the present pa- 
per will be a crucial preliminary step in view of truly 
quantum effects^ 2 - that are expected to take place in such 
miniaturized systems whenever the Joscphson charging 



energy for a single polariton becomes comparable to both 
the linewidths and the hopping energy. This regime is ex- 
pected to be entered in the next generation of samples. 

In Sec. mi we introduce the effective two mode model 
and we write the motion equations describing the time 
dynamics of the polariton field amplitudes in each of the 
two wells. The phase diagram and the different insta- 
bility regimes are analyzed in Sec. IIIII for the case of a 
continuous-wave pumping. The result of a numerical in- 
tegration of the dynamical equations of the two-mode 
model is presented in Sec. IIVI for the case of a quasi- 
continuous wave pump and in Sec. [V] for the case of a 
pulsed pump: optical bistability and self-pulsing phe- 
nomena take place in the former case, Josephson oscilla- 
tions and self-trapping in the latter one. In Sec. IVI| the 
predictions of the two-mode model are compared to full 
numerical simulations of the generalized Gross-Pitaevskii 
equation. The observability of all the predicted features 
is verified using realistic parameters for coupled polariton 
boxes. Conclusions are finally drawn in Sec. I VIII 



II. THE TWO-MODE MODEL 

A widespread description of the Josephson dynamics 
in two-well system is based on an effective two-mode 
modeli 7 -. In addition to the linear coupling J and the 
cubic nonlincarity g, a coherent pumping -Fi.,2 (t) and 
loss rates 71,2 arc to be included in order to describe 
the driven-dissipative nature of the present system. The 
equations of motion for the mode amplitudes ^1,2 (t) then 
read 2 ^: 

ih^x = (hw! - Vi + g\ fa I Vi- Jfo + F^tUi) 
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In the absence of nonlincarity g = and pumping iq^ 
0, the eigenvalues of the linear equations are 



where n,- 



E+.- = 



1 



±i^/jft(^ - w 2 ) - i( 7l - 72 )] 2 + 4J 2 : (3) 

the linear coupling J splits the unperturbed levels Tiuji^ 
into a pair of mixed eigenmodes with energies £+,-• 
For zero detuning, u>\ = u>2> an d equal loss rates 71 = 
72, the energy splitting is £L — £L|_ = 2 J, and the 
two corresponding eigenmodes are a symmetric mode 
ip+ = ip s = (ipi + - 02)/v / 2 and an antisymmetric mode 

Under a symmetric pump, F±(t) = -^(t), only the sym- 
metric mode ip+ is excited, while under an antisymmetric 
pump, Fi(t) = — i*2(t) only the antisymmetric mode 
is excited. Under a pump acting only on one unperturbed 
mode that is F\{t) = F(t), F 2 = both the eigenmodes 
are excited. In what follows, we will concentrate our at- 
tention on this last case. 

This simple linear analysis is made richer by the pres- 
ence of nonlinear terms. Actually, a cubic nonlincarity 
has two main, and strictly related effects: it introduces 
intensity-dependent shifts of the effective energy levels 
and can be responsible for dynamical instabilities. Be- 
cause of the nonlincarity, the effective eigenmodes of the 
system are no longer the symmetric and the antisymmet- 
ric ones. Therefore, to refer to the actual eigenmodes of 
the system, we prefer to adopt the notation "+" and "-" 
mode. 



\ipj\ 2 defines the stationary intensity in the 



two modes. From Eqs. ([5]), it is straightforward to show 
that the mode amplitude ip* and the pump amplitude F s 
are uniquely determined for each given pair of values of 
the pump energy Tiw and of the stationary amplitude 
in the non-pumped mode 2. By rearranging Eqs. ((5|, one 
then obtains the final equations 



F* 



= J- 



T1W2 — fiut 
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JM, (7) 



from which the F s (Huj, n-i) and tpf (fiu), n 2 ) diagrams in 
the frequency-intensity (Tko, 7i2)-plane shown in Sec. IIIII 
will be obtained. 



B. Stability of the stationary solution 

The stability of the stationary solutions found in the 
previous section can be assessed evaluating the spectrum 
of small fluctuations around the stationary solution: 



^(t)=e-*"* (^ + <%(t)), 



(8) 



By linearizing the motion equation Eqs. |T][2|) around the 
stationary solution, one obtains the following linear equa- 
tions: 



A. The stationary state 

In the whole paper, we shall restrict our attention to 
the case of a pump acting on a single mode, i.e. Fj (t) = 
8j\F(t): this scheme is in fact the most interesting for 
applications and is amenable to an almost fully analytical 
treatment. 

We start by considering the steady state of the sys- 
tem under a continuous monochromatic pump, F\{t) = 
e~ tut F s , where the amplitudes of the two modes oscillate 
at the pump frequency 

i> j {t) = e- i « t ry (4) 

Substituting this ansatz into Eqs. (HE]), the following sta- 
tionarity equations are immediately obtained: 

(/'-, - nu - i'-f) rj+gn^-Jra-j+SjiF 8 = , (5) 




+ 2g \r 3 \ 2 61>j + g {^f <ty? - JS^.j. (9) 
Substituting in Eqs. @ the time evolution 

8^ 3 (t) = e- lEt/h U 3 + e iEH ' h V*, (10) 

expressed in terms of the excitation energies E and of the 
fluctuation amplitudes Uj and Vj, the problem reduces 
to the secular equation 

M -8$ = (11) 

where we have introduced the vector 6^ = (Ui V1U2V2) 
and the matrix M has the Bogoliubov form 
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in terms of the frequencies ujj = luj — uj. 

The resulting spectrum consists of four eigenvalues E a , 
a = 1, 4, corresponding to the normal modes 8^> ' a . As 
shown by Eq. (JTDJ) , if the imaginary parts of all the 
four energies are negative Imji^Q,} < 0, the fluctuation 
is damped and the stationary solution is stable. On the 
other hand, if the imaginary part of at least one eigen- 
value is non-negative lm{E a } > 0, the solution is unsta- 
ble. 

In this latter case, two situations are possible: if 
Kc{E a } = the solution is one- mode (saddle-node) un- 
stable (1M), while, if Kc{E a } > 0, the solution is para- 
metrically unstable (P). These two situations will be dis- 
cussed in detail in what follows. 



III. CONTINUOUS EXCITATION: PHASE 
DIAGRAM AND FLUCTUATION SPECTRUM 



Under continuous monochromatic pumping, a contour 
plot of the pump amplitude F s as a function of the pump 
energy Tiuj and of the intensity n 2 in the non-pumped 
mode is readily obtained from Eqs. (07]). An example of 
such plot is shown in Fig. [T] for the symmetric u>\ = uj 2 
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FIG. 1: Energy-intensity (ftw-r^) diagram. The colour scale 
corresponds to increasing values of \F S \ in logarithmic scale 
and the thin isolines are geometrically spaced contour lines. 
The thick black line separates the stability (S) regions from in- 
stability ones: depending on their character, these are marked 
by either 1M (one-mode instability) or P (parametric insta- 
bility). The dash-dot line separates a one-mode instability 
region from a parametric instability one. The dots on the 
vertical hu = 0.45 meV line correspond to the pump pa- 
rameters used later on in Figs. [7J and [8] The system pa- 
rameters are inspired from the symmetric double-box po- 
lariton traps discussed in Sec. IVII namely J — 0.5 meV, 
71 = 72 = 7 = 0.2 meV, g = 1.1 x 10~ 3 meV, loi = uj2- 



The phase diagram is determined from the stability of 
the stationary state The black thick lines mark the 
contours between the regions of stability and the regions 
of instability, as obtained by solving the linearized prob- 
lem Eq. (p~2|) . Here, regions of one-mode instability or 
parametric instability are indicated by respectively 1M 
and P following the definitions introduced in Sec. Ill Bl 

For very small pump amplitudes, two resonances are 
clearly visible in Fig. Q] As explained in Sec. HH they 
corresponds to the ip± eigenmodes of the linearly coupled 
system and lye at exactly Tiuj — ± J. The linewidth has 
been taken smaller than the linear coupling 7/2 < J, so 
the corresponding lines are well distinct. 

For increasing values of the pump amplitude, the reso- 
nances of the system are modified by the nonlinearity (see 
the scheme in Fig. ^ . At moderate pump amplitudes for 
which gri\ 2 < J, 7, the main effect of the nonlinearity is 
a blue-shift of the two resonances. On the other hand, 
several different instability mechanisms can take place at 
larger pump amplitudes depending on the pump energy 

hui. 



A. Optical Limiter 

For Tilo < — J, the pump energy lies below the two en- 
ergy levels: this configuration, often called optical limiter 
(OL) in the literature^, is stable for all pump intensities. 
As shown in the sketch in Fig. [2ja), the effective energy 
levels are in fact pushed even further off resonance from 
the pump by the nonlinear blue-shifts. 



B. Optical bistability 

For a pump energy just above the lower energy + mode 
(i.e. Tilj <; — J), the nonlinear shift is able to push the 
effective energy level into resonance with the pump (see 
Fig. EJJb)) and give rise to a single mode (1M), saddlc- 
node^i instability. Analogous behaviour takes place for 
pump energies just above the higher energy — mode (i.e. 
for hu > J). 

The occurrence of regions of one-mode instability for 
energies higher than a mode resonance is a well stud- 
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FIG. 2: Scheme of the effective energy levels E± of a two 
mode model under a continuous-wave pump. The vertical 
dashed arrows indicate the nonlinear blue-shift. The horizon- 
tal solid arrows indicates the pump energy hu> in the different 
regimes discussed in the text: (a) Optical limiter; (b) One- 
mode instability; (c) Parametric instability; (d) High-energy 
pumping. 
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FIG. 3: Intensities m (panel (a)) and ri2 (panel (b)) as a func- 
tion of the pump amplitude. Pump energy hui = 0.45 meV. 
The arrows highlight the hysteresis cycle due to optical bista- 
bility. Same system parameters as in Fig. [1] 



ied subject in the general literature on instabilities^. 
Concerning nonlinear optical systems, it has been ex- 
tensively studied both in the simplest case of single cavi- 
ties^, as well as in more complex cases of coupled optical 
cavities^ and POs^L28. 

As a general feature^, one-mode instabilities of this 
kind often give rise to bistable behaviors, i.e. the coex- 
istence of several stable solutions for the same values of 
the pump energy and amplitude. An example of this be- 
haviour is shown in Fig. [3] where the dependence of the 
intensities nifi on the pump amplitude is plotted for a 
pump energy just above the lower resonance. 

An hysteresis cycle is apparent: as the pump ampli- 
tude F increases from zero, the system moves along the 
lower branch of stable solutions until its end point is 
reached. Only at this point the system jumps on the 
upper branch. If the pump amplitude is then decreased, 
the system keeps moving along the upper branch of sta- 
ble solution until its end point is reached, where it jumps 
back to the lower branch. 



C. Parametric instability 

For a pump energy between the two resonances, a 
parametric instability appears. When the pump energy 
equals the average of the effective energies (E + + E-)/2, 
the parametric process^ sketched in Fig. [2jc) where 
the pump field creates a signal + and an idler — fields 
becomes resonant. This happens within the window 
7j < fko < J + gn,2, where the stationary solutions are 
stable for small pump amplitudes, but become paramet- 
rically unstable as soon as the parametric gain is able 
to overcome the losses gn 2 > jj. Note that the effect 
of m can be neglected here, as in the considered energy 



window one has n\ <C n 2 (see Fig. [5J. In the dynami- 
cal systems language, such an instability is called a Hopf 
bifurcation^!. 

The strong amplification of fluctuations around the 
stationary solution eventually results in a self-pulsing dy- 
namics where the system keeps on oscillating for indefi- 
nite times. From a different point of view, these oscilla- 
tions can be seen as the result of the interference between 
three fields at different frequencies, i.e. the pump, sig- 
nal, and idler fields of the parametric oscillator^^. This 
behavior will be discussed in better detail in the next 
section. 



D. High-energy region 

For pump energy exceeding the energy of both reso- 
nances, several instability regions are expected to appear 
as a consequence of the complex interplay of single-mode 
1M and parametric P instabilities. Both effective energy 
levels eventually cross the pump energy, as well as the 
parametric resonance [Fig. |2^d)]. The diagram in Fig. [T] 
is therefore much richer in this window: for a given pump 
energy Tiuj and increasing values of n 2l three regions of 
one-mode instability and three regions of parametric in- 
stability can be identified, as well as a thin stability re- 
gion in between the first 1M and P regions. 

To understand the origin of the different regimes, the 
intensities n\. 2 are plotted in Fig.Uas a function of for 
tiuj = 1.5 meV. The behaviour is quite complex, yet can 
be analytically interpreted from Eq. ([6]), which indeed 
gives 



Til = J 



(gn 2 - JvjS) + -r 



n 2 



(13) 
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FIG. 4: Intensity m as a function of ri2 for a high pump energy 
hu — 1.5 meV (dashed line). The solid line corresponds to 
ri2 and is a guide for the eye in order to identify the regions 
where ni > «2, «i 77,2 or ni 3> n\. Separation between 
instability regions of different (1M or P) kinds are indicated 
by the vertical dotted lines. Same system parameters as in 

Fig.m 



For a pump energy lying between the two blue-shifted 
resonances, i.e. —J + gn 2 < Tluj < J + gn 2 , the intensity 
n 2 of the non-pumped mode is larger than the intensity 
m of the pumped one, i.e. n 2 > n±. Conversely, for 
either lower (%u> < — J + gn 2 ) or larger (Hlu > J + gn 2 ) 
pump energies, the pumped mode has a larger intensity. 
This trend is illustrated in the ri\ jn 2 plot shown in Fig. [5] 
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FIG. 5: Logarithmic gray-scale plot of the ratio 711/712 as 
a function of the pump energy and of the value ri2. Same 
parameters as in Fig. [T] 

First we investigate the region around the stability 
tongue extending at hoj > 1.5 meV for relatively small 
values of F. While the energy of the pumped 1 mode is 
significantly blue-shifted by the nonlinear term, the non- 
pumped 2 mode remains almost empty (see Fig. 2]). The 
two 1,2 modes are then far in energy, so the effect of 
the linear coupling is strongly suppressed. The nearby 
1M instability region then corresponds to the bistability 
loop for a pump close to resonance with the higher en- 
ergy mode, which in this region basically coincides with 
the blue-shifted 1 mode. The physics is analogous for the 
second 1M instability region located just above. 

The third 1M instability at much larger pump ampli- 



tudes corresponds to the opposite situation where the 
2 mode has been shifted above Tiuj and consequently 
n 2 ^ ni according to (fT5|) . The bistability loop then 
involves the lower resonance, which in this regime basi- 
cally corresponds to the unperturbed 1 mode. 

We finally consider the intervals of parametric insta- 
bility. The first and the third intervals correspond to 
resonant scattering processes taking place in a regime 
where the modes are effectively decoupled as n\ 3> n 2 
or n 2 3> ni, respectively. In these regimes, the two ± 
cigenmodes essentially coincide with the 1, 2 modes. The 
second interval corresponds instead to an intermediate 
regime where n\ ss n 2 , and the two ± eigenmodes are a 
superposition of both 1, 2 modes. 



E. The spectrum of fluctuations around the 
stationary solution 

The stability properties of the stationary solution dis- 
cussed in the previous section are further illustrated by 
looking at the eigenvalues of the Bogoliubov linearized 
theory (fT2")) of small fluctuations around the stationary 
stated. These are plotted in Fig. [5] as a function of the 
intensity n 2 of the non-pumped mode for the case of a 
pump energy hu chosen between the ± eigenmodes of the 
unloaded system. 
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FIG. 6: Real (a) and imaginary (b) part of the linearized 
spectrum around the steady state solution as a function of 
712. Stable regions correspond to the imaginary part of the 
spectrum being negative. Pump frequency hco = 0.45 meV. 
The vertical lines indicate the 712 values used in the panels 
(a-d) of Figs. [7] and [5] Same system parameters as in Fig. [T] 

In the linear 711,2 — > regime, the frequencies and 
damping rates tend to the ones Rc{i?3,i} = — Reji^} = 
±J/2 — hu;, — Im{i?i,2,3,4} = 7/2 of the unloaded sys- 
tem. As a consequence of interactions, the frequency 
Kc{Ei 3 } of the positive-weighted Bogoliubov modes are 
blue-shifted for growing intensities, while the ones E 2 4 
of the negative-weighted ones are red-shifted: this makes 
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them to pairwise intersect at some value of the inten- 
sity. Here, each pair collapses onto values Re{-E 2 (4)} = 
Rc{_E!(3)} opposite in sign. Correspondingly, the damp- 
ing rate increases —lm{Ei} = — ImjEi} > 7/2 for two 
of them, while it decreases — Im{i?2} = — Im{£ , 3} < 7/2 
for the two others, possibly giving rise to a dynamical in- 
stability, as in the case displayed in the figure. The fact 
that frequencies of the modes involved in the instability 
are non-zero is a signature of the parametric nature of 
the instability. 

For larger intensities, the frequencies split again within 
a narrow intensity interval where the damping rate goes 
back to 7/2, but another instability region occurs at even 
higher intensities as a consequence of the intersection 
Rc{E\} = Rc{Ei}: while the the imaginary parts of the 
2,3 modes stay unchanged, the ones of the 1,4 modes 
are split and dynamical instability is signalled by one of 
them becoming positive. Since the unstable mode has 
zero frequency, the instability has the one-mode charac- 
ter typical of optical bistability loops. 

For even larger intensities, the stationary state be- 
comes stable again: because of the large blue- (red-) shift 
of the positive (negative)-weightcd Bogoliubov modes, no 
further intersections of the mode frequencies can in fact 
occur. 



IV. CONTINUOUS PUMP: BISTABILITY, 
SELF-PULSING, AND RESPONSE TO A PROBE 

The stationary states and the stability regions iden- 
tified in the previous section are a good starting point 
for the dynamical study of the system that we carry out 
in the present section by numerically solving Eqs. (Q} 
[2]). We first investigate the onset of the steady state 
when the pump intensity is slowly increased in time to 
its asymptotic value. Then we characterize the response 
of the system in its steady state to an additional probe: 
this provides a simple and effective way of measuring the 
frequencies and damping rates of the Bogoliubov modes 
in the different regimes. 

The quasi-continuous pump is assumed to have a 
smooth temporal profile of the form 



enough for the pulse to encompass all the relevant spec- 
tral features, i.e. all the 4 Bogoliubov modes shown in 
Fig. El 
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t(ps) 

FIG. 7: Time-evolution of the n\(t) (dashed line) and ri2(t) 
(solid line) intensities for growing pump amplitudes F max = 
lmeV (a), 5meV (b) 8meV (c), 20meV (d). Pump energy 
hw = 0.45 meV. Pump switch-on time r = 100 ps. Probe 
amplitude f° = 10~ 2 Fmax <C F ma x, probe duration r g = 
0.3 ps <C irh/J, probe delay to = 400 ps 3> t. Same system 
parameters as in Fig. [1] The pump parameters used for panels 
(a-d) correspond to the red dots in Fig. [T] 
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For very long switch-on time, the system evolves in a 
quasi-static way through a sequence of stable stationary 
states. 

Once the system has got to its asymptotic stationary 
solution, a weak and short probe pulse is applied onto 
mode 1. Its temporal shape is assumed to be a Gaussian 



(15) 



its central frequency coincides with the one of the con- 
tinuous pump, and its duration r ff is chosen to be short 



Numerical predictions for ni^it) are shown in Fig. [7] 
for the different regimes. The quasi-cw pump energy Tiuj 
is taken to be between the linear resonance peaks and in- 
creasing values of the pump amplitude F are chosen for 
the different panels (see dots in Fig. [I} . The correspond- 
ing spectra shown in Fig.[H]are obtained by Fourier trans- 
form of ipi^it)- In order to eliminate complications due 
to the switch-on dynamics, we have restricted the Fourier 
transform to the temporal window following the arrival 
of the probe pulse. The central (5-peak corresponds to 
the pump frequency. 
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FIG. 8: Fourier transform spectra of tfii (dashed line) 
and tp2 (solid line) for growing pump amplitudes F max = 
1, 5, 8, 20 meV. The central 5-peak is at the pump frequency. 
The arrows indicate the frequencies of the Bogoliubov modes 
with the same colour code as in Fig. [6] Parameters of the 
(a-d) panels correspond to the ones in Fig. [7] 



A. Stable regime 

For small values of the asymptotic pump amplitude 
Fmax [Fig.EIa) and (b)], the evolution during the switch- 
on time smoothly leads the system to the asymptotic 
stationary state. As this solution is stable, the response 
to the probe pulse gets quickly damped within a time- 
scale of the order of ten picosecond. While for very small 
intensities (Fig. |7fa)) the response consists of damped 
oscillations at a single frequency, for larger intensities 
(Fig. E]Jb)) relaxation is more complex and involves in- 
terference of more frequencies. 

This difference is apparent in the corresponding spec- 
tra shown in Fig. [Hfa,b) which are to be compared to 
the Bogoliubov modes shown in Fig. [51 Well inside the 
stability region, only the positive-weighted Bogoliubov 
modes with a significant U component [see Eq. (|10[) ] are 
in fact visible. On the other hand, when the parametric 
instability region is approached, the normal components 
U are significant for all modes and all the four frequen- 



cies become then visible in the spectrum (see Fig.[8]Jb)). 
As usual, the finite linewidth of the peaks is fixed by 
the finite and negative imaginary part of the Bogoliubov 
modes, i.e. by their damping rate. In the present case, 
this is the same for all Bogoliubov modes. 



B. Parametric instability 

The physics is richer in Figs. \T[c) and [5Jc) where a 
larger pump amplitude F max = 8 meV is considered: in 
this case, the asymptotic stationary state is in fact para- 
metrically unstable. As the pump intensity is increased 
very smoothly, the system adiabatically follows the sta- 
tionary state at the instantaneous value of F even in the 
unstable region. However, the arrival of the probe pulse 
speeds up the onset of the instability: the perturbation 
that it induces is quickly amplified until the system gets 
to the self-pulsing regime where undamped periodic oscil- 
lations take place for indefinite time. Their frequency is 
close to the one of the linear Bogoliubov mode getting un- 
stable. Correspondingly, two 5-peaks appear in the spec- 
trum shown Fig. [SJc) at energies E — fiu — ±0.52 meV. 
The weaker 5-peaks at harmonic frequencies contribute 
to the quite complex waveform of the self-pulsing oscil- 
lations in time shown in Fig. [7Jc) ■ 



C. One- mode instability 

To highlight hysteresis phenomena, we now choose an 
asymptotic value of the pump amplitude above the jump- 
up threshold of the bistability loop shown in Fig. [3l The 
switch-on then does not take place in a smooth way and 
the system has to perform a sudden jump when the end- 
point of the lower branch is reached: among the many 
complex behaviour that may take placed, the behaviour 
of the present system is the simplest: it jumps to the 
higher-intensity branch of the bistable loop. Once these 
upper branch is reached, the system is stable again and 
the response to the probe pulse is qualitatively similar 
to the case of panels (a,b). The only difference is the 
higher frequency of the oscillations, and the presence of 
three different excitation frequencies which contribute to 
the complex relaxation dynamics: one negative-weighted 
Bogoliubov mode has in fact a significant U component. 



V. PULSED EXCITATION: JOSEPHSON 
OSCILLATIONS AND SELF-TRAPPING 

After having investigated the behaviour of the system 
under a continuous pump, it is now interesting to look at 
the case where only a pulsed pump is applied to the sys- 
tem. Specifically, we numerically solve Eqs. using 
a Gaussian temporal profile for the pump 



F{t*) — Fmax & 



<t-tpf/< 



(16) 
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The duration t p of the pump pulse is taken to be very 
short as compared to all time scales of the system dy- 
namics: the system is then almost istantaneously excited 
by a sudden kick, and then let evolve and relax with- 
out any further pumping. The results are summarized in 
Fig. [9ja,c,e,g) in the time domain, while the correspond- 
ing Fourier spectra are shown in Fig. [5Jb,d,f,h). 




10 20 -2 o 2 



t (ps) E (meV) 

FIG. 9: Intensity dynamics (a,c,e,g) and corresponding 
Fourier spectra (b,d,f,h) under a Gaussian pulsed pump (Eq. 
(|16l> ) at time t p = 3ps of duration r p = 0.2 ps. Peak pump 
amplitude F ma x = 1 meV (a,b), 50 meV (c,d), 75 meV (e,f) 
and 100 meV (g,h). Blue (red) lines corresponds to the 1 (2) 
modes. Spectra are obtained by Fourier transforming the sig- 
nals in the whole time interval. Same system parameters as 
in Fig. [JJ 

For low pump amplitudes in the linear regime, the 
intensities in the two wells manifest Josephson-like os- 
cillations [Fig. [5(a)]. The pump pulse creates in fact a 
localized excitation in the 1 mode, which is a superpo- 
sition of the symmetric and anti-symmetric eigenmodes 
of the system. Because of their energy splitting, the sys- 
tem shows complete Josephson oscillations with a period 
nh/ J, which then damp out at a rate 7 under the effect 
of losses. Correspondingly, the Fourier spectrum shown 
in Fig. [FJb) is characterized by a pair of resonance peaks 
split by 2 J. 

For stronger pump amplitudes, the instantaneous in- 
tensity in the system increases, and eventually results 
in significant nonlinear effects. While the time evolu- 
tion of the intensities [Fig. [DJc)-(e)] does not appear to 
be significantly modified, nonlinear effects are visible in 
the spectra of Fig. [9jd)-(f) already at moderate intensi- 
ties as a global blue-shift of the spectrum and a signifi- 
cant increase of the width of the peaks. In particular, in 
Fig- Elf) > note how the blue-shift of the spectrum relative 



to the field ipi is larger than the blue-shift of the spectrum 
relative to ip2'- this is due to the fact that the intensity 
n\ is at short times much larger than the intensity . In 
the Fourier spectrum of ip\ , we also clearly recognize two 
peaks at J = ±0.5 meV corresponding to the eigenfrc- 
quencies of the linear dynamics that is recovered at long 
times once the intensities have dropped to small values. 
Similar peaks also contribute to the Fourier spectrum of 
i/>2, but are hardly visible in the figure, their weak inten- 
sity being hidden by the tails of the main peaks. 

The appearance of peaks at frequencies characteristic 
of the nonlinear regime is a precursor of the self-trapping 
regime that appears for stronger pump amplitudes: in 
this case, the nonlinear effects are in fact dominant in 
determining both the spectrum and the time evolution 
of the intensities [Fig. [FJg)-(h)]. 

In the early stages of the evolution when the intensity 
is the largest, nonlinear effects dramatically suppress the 
amplitude of Josephson oscillation: most of the intensity 
is in fact self-trapped in the mode 1 and the intensity of 
mode 2 oscillates around a much smaller value. As time 
goes on, the total intensity slowly drops under the effect 
of losses and eventually complete Josephson oscillations 
are recovered: the transition between the self-trapping 
regime and Josephson oscillations can be located in the 
vicinity of the time when the total intensity equals the 
critical density n\ + 712 = n c tot = 4:J/g of equilibrium 
Josephson systems^!. The presence of losses is only re- 
sponsible for a small shift of the critical point. 

Consequences of this physics can be observed also in 
Fig. E(h) : both spectra show in fact two broad peaks cen- 
tered at high energies (i.e. at E ~ and E ~ 1.5 meV), 
which represent the two resonances of the system in 
the self-trapping regime, and two lower energy peaks at 
E = ±J = ±0.5 meV, representing the frequencies of the 
Josephson oscillations. These two latter peaks are asym- 
metric [differently from the pure linear regime displayed 
in Fig. Mb)] because the dynamics has been modified by 
the occurrence of the self-trapping regime. 



VI. MICROCAVITY POLARITON BOXES 

In this last section we show how the parameters of the 
two-mode model can be evaluated from the microscopic 
structure of a specific physical system. On one hand we 
demonstrate that all the physics discussed in the previous 
sections can actually be observed in realistic systems, on 
the other hand we confirm the quantitative validity of the 
predictions of the two-mode model by comparing them 
to a full numerical integration of the generalized Gross- 
Pitaevskii equation for polaritons 3 ^. 

Although many other configurations based on e.g. cou- 
pled DBR cavities are available to study Josephson-likc 
effects in an optical contex t 19 ! 20 ' 29 ' 30 , our attention will 
be concentrated on the specific case of double-well po- 
lariton traps. Such a system was recently realize d 12 ' 13 
and combines the strong nonlincarity due to the exci- 
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tonic component of the polariton to the possibility of a 
micron-scale spatial confinement by laterally patterning 
the thickness of the cavity layer. From the point of view 
of Josephson physics, this geometry is very attractive as 
it allows independent collection of light emitted from the 
two boxes and preserves the signal from being covered by 
the incident laser field. 



From the non-equilibrium GPE to the two 
mode model 



The dynamics of the macroscopic polariton field $(r, t) 
is described at mean-field level by a non-equilibrium gen- 
eralization of the Gross-Pitaevskii equation (GPE) of the 



2m n 



U ext (r)-i^ )¥(r,i) 



+ v \*(r,t)\ 2 V(r,t)+f(r,t) . (17) 



Trip is the effective mass of the lower polariton, 7 is the 
decay rate, U ex t is the trapping potential, v is the effec- 
tive polariton mutual interactio n 31 ' 32 ! 33 and /(r, t) is the 
amplitude of the coherent pump field. 

In this paper, we consider the case of a trapping poten- 
tial Uexti?) formed by two adjacent wells, as displayed in 
Fig.rUa): 

U ext = -U Q 9 (L y - \y\) [9 (L x + 5/2 - x) 9 (x - 5/2) 
+ 9 (L x + 5/2 + x) 6 (-x - 5/2)] . (18) 

For this potential, the fundamental mode of energy E gs is 
described by an cigenfunction 4> gs (r) which is symmetric 
in the two wells (see Fig. rXUT b)). while the first excited 
mode of energy E exc corresponds to an antisymmetric 
eigenfunction (f> eX c(r) (see Fig. [TUTc)). Without loss of 
generality, both these functions can be taken real. 

The pump field is considered to be monochromatic, 
with a gaussian spatial profile 



f(r,t) =27ra 2 /oe- 



(19) 



centered on the first well at ri = ((L x +8)/2,L y /2) t . Pro- 
vided the pump energy is close to E gs and E exc and all 
other excited modes of the trapping potential are at much 
higher energies, the coherent field ^(r, t) can be safely 
written in the two-mode limitii as a time-dependent su- 
perposition of the two lower energy modes ffs (r) and 
4> exc (r) only. 

For the present case, it is useful to write the superpo- 
sition in the form 



* (r, t) = V>i (f) 0i (r) + V> 2 (t) 02 (r) 

where the Wannier-like functions 
1 



0i,2 W 



71 



(r) ± cj) exc (r) 



(20) 



(21) 




FIG. 10: Trapping potential U ex t (Eq. (Tl8)) due to two ad- 
jacent rectangular wells (a), and the fundamental ss (r) (b) 
and first excited <f> eX c{r) (c) modes of the corresponding equi- 
librium GPE (|3ip . The spatial profiles are displayed in grey 
tones. 



are mostly localized in each well and orthogonal to each 
other. 

By substituting Eq. ([ST))) in Eq. (jTTJ) and then project- 
ing onto the lowest states, two coupled dynamical equa- 
tions for the amplitudes f/TsW are obtained of the form: 

ihtpj = (huij — £7/2) tpj + g ipj — J1P3-J 

+ 6[(2|^| 2 + |V3- J -| 2 )V'3-i + ^-^ 2 ] 

+ c + r^l-j) + Fj(t). (22) 

Linear dynamics is summarized by the diagonal term 

(23) 



llUj — 7. {Egs + E exc ) 



and the linear hopping coefficient 



J = - (E eXC - Eg S ) . 

Pumping is described by 



(24) 



(25) 



while nonlinear effects are described by the three coupling 



F j {t)= / dr^(r)/(r,t), 
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coefficients 



9 
b 



rfr^i, 2 (r)] 4 
dr0 2 ,i(r) [0i,2(r)] c 



c = « / dr[0i(r)] 2 [^(r)] : 



(26) 
(27) 
(28) 



For typical geometries and for moderate intensities, the 
coefficients b and c are much smaller than the other quan- 
tities and can be safely neglected. Within this approx- 
imation, Eq. (f2"2"|) reduces to the two mode model Eqs. 
([I][2]) used in the previous Sections. 



B. Comparison with the 2-mode model 

In order to verify the validity of the two-mode approx- 
imation, a numerical integration of the full GPE Eq. JT7|) 
can be performed and then quantitatively compared to 
the predictions of the two-mode model. 

For this comparison, realistic parameters for typical 
polariton boxes in GaAs based microcavitie o 12 ! 13 are 
used, that is a trapping potential depth Uq — 5 meV, 
lateral box sizes L x = 1 /xm and L y = 2 /im, a polariton 
mass m p = 7 x 10 mo (too being the electron mass) and 
a nonlinear coupling constant v = 2 x 10~ 3 meV/im 2 . We 
assume the separation 8 = 0.5 /im between the two wells. 
The corresponding potential U ex t and the functions rel- 
ative to the first two GP modes are displayed in Fig. [TO] 
A spatial width a = 0.5 /im is taken for the pump spot. 

The time-evolution of the polariton occupation 
ni,2(t) = | "01,2 | 2 in the two wells is plotted in Fig. QT] 
In the present GPE framework, the occupations ni.2 of 
each of the two wells are defined by the spatial integrals 



(29) 
(30) 



rn{t) = J dxjdy |*(r,t)| 2 
n 2 (t) = f dx fdy |f(r,t)| 2 . 

J — oo J 



Two different parameter choices are made in the two pan- 
els of Fig. QTJ In (a) , the system tends to a stable station- 
ary state, while self-pulsing oscillations are visible in (b). 
The qualitative agreement with respectively panels (a,b) 
and (c) of Fig. [7] is apparent. Note that the switch-on 
time of the pump considered in Fig.[Tl]is not long enough 
to guarantee a quasi-static evolution of the system even 
in the unstable region: differently from Fig.[7Jc), the self- 
pulsing oscillations arc immediately visible without the 
need of a perturbation seed. 

As the linear coupling J is very sensitive to the shape 
of the wavefunctions in the barrier, interactions may af- 
fect it quite significantly, spoiling the quantitative agree- 
ment with the two-mode model. In order to get a good 
quantitative agreement in all regimes, a better approxi- 
mation can be adopted for the localized wavefunctions us- 
ing for 4> gs (r) and <p exc (r) the two lowest-energy solutions 
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FIG. 11: Dynamical evolution of m and ni (ni > m), as ob- 
tained from the full numerical calculation of GPE (thin lines) 
and from the two-mode model (thick lines). Panel (a): GPE 
simulations with pump energy hu = —2.8 meV and pump 
amplitude fo = 5 meV; parameters of the corresponding ef- 
fective two-mode model: ?iwi,2 = —3.15 meV, J = 0.45 meV, 
g = 10~ 3 meV and Ft = 5.8 meV. Panel (b): GPE simula- 
tions with pump energy hui = — 2.8 meV and pump amplitude 
fo = 7.5 meV; parameters of the corresponding effective two- 



mode model: Hujj = 
and Fi = 8.5 meV. 



-3.1 meV, J = 0.5 meV, g = meV 



of the time-independent, equilibrium Gross-Pitaevskii 
equation: 



+ U ext {r) + vN |0(r)| 2 ) ^(r)=#^(r). 



ft 2 V 2 

2m p 

(31) 

The total polariton number N = n\ + n-j has to be ob- 
tained from the long-time limit t — > oo of the full GP 
equation once the solution has come to their asymptotic 
steady state, or by averaging over the period of the self- 
pulsing oscillations. Correspondingly, Eqs. (|23]) and (|24|) 
have to be substituted by the N-dependent equations 



~ J dr([0 ss (r)] 4 + [0 e;EC (r)] 4 ) (32) 



and 



(33) 



respectively. As expected, the main effect of interactions 
is to slightly lift the bottom of the two wells, so to effec- 
tively reduce the barrier height and enhance tunneling. 

The result of such a procedure is also shown in Fig.fTTI 
the overall qualitative agreement is good. From a quanti- 
tative point of view, the agreement is always excellent in 
the stable regime of panel (a) , while some discrepancies 
arc visible in panel (b), in particular at short times be- 
fore the self-pulsing sets in. This can be expected, as the 
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parameters of the two-mode model extracted from the 
late time dynamics slightly differ from what one would 
get from the early stages. Furthermore, the spatial pro- 
file of the wavefunction has a significant variation in time 
during the self-pulsing dynamics. 

VII. CONCLUSIONS 

In this paper we have studied the two-mode dynamics 
of spatially coupled polariton boxes under a coherent ex- 
ternal pumping and we have shown that this provides an 
interesting non-equilibrium optical generalization of the 
well-known Josephson effect of weakly coupled superflu- 
ids and superconductors. 

For a continuous wave pumping, a phase diagram has 
been obtained which summarizes the steady state of the 
system as a function of pump intensity and frequency. 
Stable and unstable regions have been identified; one- 



mode and parametric instabilities have been shown to 
be intrinsccally related to respectively optical bistability 
and self-pulsing effects. The response of the system to 
an additional probe provides unambiguous information 
on the Bogoliubov modes around the stationary state. 

For a short pump pulse, the crossover from a Joseph- 
son oscillations regime to a self-trapping one has been 
characterized as a function of the pump intensity. Pe- 
culiar features due to the non-equilibrium nature of the 
system have been pointed out. 

The validity of the two-mode model and the actual 
observability of the predicted effects has been verified on 
the basis of the non-equilibrium Gross-Pitaevskii equa- 
tion for polaritons in a double-well trap potential using 
parameters inspired by recent experiments. 
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